Interaction Effects in Topological Superconducting Wires Supporting Majorana 

Fermions 



E. M. Stoudenmire, 1 Jason Alicea, 1 Oleg A. Starykh, 2 and Matthew P. A. Fisher 3 

1 Department of Physics and Astronomy, University of California, Irvine, CA 92697 
2 Department of Physics and Astronomy, University of Utah, Salt Lake City, Utah 84112 
3 Department of Physics, University of California, Santa Barbara, CA 93106 
(Dated: January 20, 2013) 

Among the broad spectrum of systems predicted to exhibit topological superconductivity and 
Majorana fermions, one-dimensional wires with strong spin-orbit coupling provide one of the most 
promising experimental candidates. Here we investigate the fate of the topological superconducting 
phase in such wires when repulsive interactions are present. Using a combination of Density Matrix 
Renormalization Group, bosonization, and Hartree-Fock techniques, we demonstrate that while 
interactions degrade the bulk gap — consistent with recent results of Gangadharaiah et al. — they 
also greatly expand the parameter range over which the topological phase arises. In particular, we 
show that with interactions this phase can be accessed over a broader chemical potential window, 
thereby leading to greater immunity against disorder-induced chemical potential fluctuations in the 
wire. We also suggest that in certain wires strong interactions may allow Majorana fermions to be 
generated without requiring a magnetic field. 

PACS numbers: 73.21.Hb, 71.10.Pm, 74.78.Fk 



I. INTRODUCTION 

The search for Majorana fermions in condensed matter 
systems ispresently being pursued with great vigor and 
optimism driven to a large extent by the surge in real- 
istic proposals^"22l f or stabilizing topological phases sup- 
porting these exotic particles. Among the numerous plat- 
forms proposed, one-dimensional wires with strong spin- 
orbit interactions have e merg ed as a top contender for 
hosting Majorana modes P^*' These systems possess the 
virtue that one can generate a topological superconduct- 
ing state in a relatively simple manner, using only an ap- 
plied magnetic field and proximity-induced Cooper pair- 
ing (which has already been achieved^) . Despite the fact 
that exchange statistics is ill-defined in one dimension, 
forming wire networks allows one to harness th e non- 
Abelian braiding statistics of Majorana fermionj 5 * 24 * 25 ! 
and hence their potential for quantum computationpSHUl 
even in this setting. Thus the discovery of Majorana 
fermions in one-dimensional wires could ultimately lead 
to profound technological innovations. 

Much of the prior work on topological superconductiv- 
ity in one-dimensional wires has focused on clean systems 
in the non- interacting limit. A thorough understand- 
ing of the conditions required to experimentally observe 
and manipulate Majorana modes in these systems will 
likely require detailed studies incorporating both disor- 
der and Coulomb repulsion, which poses a challenging 
problem. Several recent works have considered the effects 
of disorder on the topological phase in a non-interacting 
wirej 30 " 38 ! interactions, however, have received less at- 
tention, even in the clean limit. There is a very nat- 
ural reason to expect that Coulomb effects may influ- 
ence the topological phase and its associated Majorana 
modes in important quantitative and perhaps even quali- 
tative ways. Namely, in the absence of proximity- induced 



superconductivity, the one-dimensional wires of inter- 
est in fact constitute repulsively interacting Luttinger 
liquids whose properties can d epart dramatically from 
those of a non-interacting wire.^EHl The interplay be- 
tween Coulomb interactions, magnetism, and proximity- 
induced pairing in such a Luttinger liquid is far from 
obvious a priori and very important to clarify for exper- 
iment. 

Gangadharaiah et al. recently made an important ad- 
vance in this direction using a field-theoretic approach.^ 
More precisely, these authors addressed interaction ef- 
fects in the limit where the Zeeman field is sufficiently 
large that only one of the spin-orbit-split bands plays a 
role. Here the problem maps onto 'spinless' fermions and 
can be studied analytically using bosonization. Repulsive 
interactions were found to suppress the pairing-induced 
bulk gap in the wire, even destroying the gap altogether 
beyond a critical interaction strength^ — a rather phys- 
ical but unfortunate finding from an experimental view- 
point. Very recent work by Sela et al. reached similar 
conclusions in a model relevant for quantum spin Hall 
edges I SI 

Here we follow up on the study of Gangadharaiah et 
al, incorporating Coulomb repulsion at both weak and 
strong Zeeman fields using a combination of Density Ma- 
trix Renormalization Group (DMRG) , Hartree-Fock, and 
bosonization methods. Our DMRG simulations indi- 
cate that while interactions indeed suppress the bulk gap 
throughout, they also lead to a second, experimentally 
beneficial effect — topological superconductivity and Ma- 
jorana fermions can be accessed at weaker magnetic fields 
and over a broader chemical potential window compared 
to the non- interacting limit. Although we do not explic- 
itly include disorder here, this result strongly suggests 
that interactions lead to greater immunity against chem- 
ical potential fluctuations in the wire, which is highly 
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desirable given that this type of disorder is likely to pro- 
vide one of the main obstacles to observing Majorana 
fermions in this setting. 

With moderately strong interactions (compared to the 
Zeeman field and pairing energies) , we show that Hartree- 
Fock theory in fact performs remarkably well, capturing 
the DMRG phase diagram semi-quantitatively. At strong 
interactions where Hartree-Fock proves inadequate, we 
employ bosonization to address the fate of the topolog- 
ical phase in a limit complementary to Ref. 41 as we 
treat both pairing and the Zeeman field as perturbations. 
Here too we show that the topological phase expands in 
phase space due to interactions (provided they are not 
sufficiently strong that pairing becomes irrelevant), and 
demonstrate consistency with DMRG simulations in this 
limit. Finally, we suggest that a wire with strong Rashba 
and Drcsselhaus couplings may be driven into a topolog- 
ical phase by interactions, even without an applied mag- 
netic field. By performing DMRG simulations on a time- 
reversal invariant toy model with strong interactions, we 
show that the formation of Majorana fermions is indeed 
possible at zero magnetic field. 

In the following section we describe the model system 
studied throughout, and set the stage by reviewing the 
physics in the non- interacting limit. Interaction effects 
are discussed in Sec. III. In Sec. IIIA we introduce the 
diagnostics we use for identifying the topological phase 
in DMRG, and present the phase diagram determined 
at finite Zeeman and pairing fields. Our Hartree-Fock 
and bosonization studies are described in Sees. IIIB and 
C, respectively. Section HID compares our bosonization 
predictions with DMRG simulations, while Sec. HIE dis- 
cusses the possibility of accessing topological supercon- 
ductivity in a time-reversal invariant system. Finally, we 
summarize our results and discuss future directions in 
Sec. IV. 



II. MODEL SYSTEM 

The system we consider is a one-dimensional semicon- 
ducting wire with Rashba spin-orbit coupling of strength 
a and an orthogonal Zeeman field V z . Electrons in the 
wire additionally inherit an s-wave pairing field A (which 
we take to be real throughout) via the proximity effect 
with a neighboring bulk superconductor. Including a lo- 
cal repulsive interaction of strength U, we employ the 
following minimal continuum Hamiltonian for the sys- 
tem: 



Ho = 



(1) 
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FIG. 1: Single-particle energy spectrum of the Hamiltonian 
"Ho when A — 0, a = 0.3, and m = 1. The dashed lines show 
the spin-orbit-split bands in the time-reversal invariant limit 
with V z = 0. Introducing a non-zero V z opens a gap in the 
spectrum at k — as depicted by the solid lines. This opens 
up a chemical potential window where the system exhibits 
only a single pair of Fermi points and thus resembles a 'spin- 
less' metal. Turning on a weak A then effectively p-wave pairs 
the 'spinless metal', generating a topological superconducting 
state supporting Majorana fermions. 



Here x is the coordinate along the wire, m is the effective 
mass for the electrons, and a x,y,z are Pauli matrices. The 
operator ip\ creates an electron with spin s =f, J,, while 
p s = tpli/jg is the density operator. Throughout, spin 
indices are summed whenever suppressed. 

Before turning to interaction effects, it will prove very 
useful to revie w the phase diagram of % in the non- 
interacting caseP^j3 Consider first the time-reversal- 
invariant limit with V z = A = 0. Here a favors align- 
ing the electron spins along or against the y-direction, 
depending on the momentum. This produces the band 
structure shown by the dashed lines in Fig. [lj where the 
spins orient along the +y and —y directions in the left 
and right parabolas, respectively. The salient feature of 
these bands is the generic presence of four Fermi points in 
the spectrum, for arbitrary values of p that are above the 
band minimum. Turning onA/0 pairs k and —k states 
at both sets of Fermi points, producing an ordinary, non- 
topological superconductor. The system always has a 
unique ground state here which is separated by a finite 
energy from all excited states. Note that as long as time- 
reversal symmetry remains unbroken, spin-orbit coupling 
plays no qualitative role, and the Hamiltonian can be 
smoothly deformed into that of a conventional ID su- 
perconductor with a = (albeit with long-range order 
induced by the proximity effect). 

Restoring V z ^ breaks time-reversal symmetry and 
cants the spins along the z direction. This modifies the 
band structure in a crucial way; as shown by the solid 
lines in Fig. [TJ a gap in the spectrum opens at k = 0. In 
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this case there exists a chemical potential window where 
the wire exhibits only two Fermi points. Here it is con- 
ceptually very useful to imagine projecting out the upper 
unoccupied band, and then viewing the wire as a 'spin- 
less' metal. Turning on a weak (compared to V z ) s-wave 
pair field A induces p-wave pairing between the fermions 
of the 'spinless' metal, driving the system into a topolog- 
ically non-trivial superconductor. In a finite system with 
open boundary conditions, exponentially localized Majo- 
rana modes 71 and 72 appear at the left and right ends of 
the wire. Up to exponentially small corrections in their 
separation, 71 and 72 represent zero-energy modes and 
therefore lead to a two-fold ground state degeneracy, in 
contrast to the situation for the ordinary superconductor 
discussed above. More precisely, if |0) is the ground state 
with even fermion parity, then |1) oc 7j|0) is the ground 
state with odd parity. All other states are separated by 
a finite energy gap from these two. 

This topological phase in fact connects smoothly to 
Kitaev's exactly solvable model for a p-wave supercon- 
ducting chain, where the appearance of Majorana modes 
is very transparent.^ Note that spin-orbit coupling is ab- 
solutely crucial for the onset of this phase in the wire. In 
particular, when a^O the spins at the Fermi momenta 
fci? and — &p in the lower band are not parallel, which 
allows the spin-singlet pair field A to have a nontrivial 
effect on the 'spinless' metal. It will be useful to keep in 
mind, however, that as V z /a increases the spins become 
increasingly aligned and hence A becomes correspond- 
ingly less effective at gapping the 'spinless' metal. 



resulting lattice model reads 



A quantitative analysis^^D f the Bogoliubov-de 
Gennes equation with [7 = reveals that the topolog- 
ical superconducting state exists only when the following 
criteria are satisfied: 



Vz > A 



-^V z 2 -A 2 <A1 < ^/V 2 -A 2 



(4) 



Outside of this regime, the system can not be viewed as 
a 'spinless' metal either because the chemical potential 
lies outside of the Zeeman-induced gap or because A is 
sufficiently large that the upper and lower bands strongly 
couple to one another. In either case the ordinary super- 
conducting state emerges. 

One of the primary goals of this paper is understanding 
how this picture evolves when (possibly strong) interac- 
tions are present. We will attack the problem with a 
combination of DMRG, Hartree-Fock, and bosonization 
methods. To simulate the wire model using DMRG, we 
construct a lattice model that maps onto the continuum 
Hamiltonian above in the low-density limit. This can be 
done in the usual way by writing T-L in momentum space 
and replacing k 2 — > 2(1 — cosfc), k — > sinfc, J k — > J2 fc , 



H 

Ho 



Ho + Hj 
E 



(5) 



\{c]cj+i + h.c.) - (p-t)c\c 



-(icja v c j+1 + h.c.) + V z c]a z c 3 



+ M c rt c ji + h.c.) 



H, 
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(6) 
(7) 



where 



- J 



Note that upon displaying h 



and the lattice constant a, the hopping strength is 
t — h 2 1 (ma ) while the continuum and lattice definitions 
of the spin-orbit coupling constant are related through 
aiatt = \ acont- After setting h = a = 1 in what follows, 
t = 1/m and «i a tt = «cont so we employ the same sym- 
bol a to denote the spin-orbit strength in both models. 
We will distinguish these for clarity only when discussing 
energy scales in the Hamiltonians momentarily. When 
U = 0, the topological superconducting state again re- 
quires the criteria in Eq. Q , exactly as in the continuum 
model. 

Since the Hamiltonians above contain many pa- 
rameters, it is worth mentioning briefly the rough 
energy scales that are expected to be relevant in the 
problem. Semiconductor effective masses are typically 
around m ~ 0.05m e , where m e is the bare electron 
mass. Assuming the Rashba spin-orbit coupling in the 
wire is comparable to the values in two-dimensional 
quantum wells featuring heavy elements, we expect 
acont ~ 10 4 — 10 5 m/s. These estimates yield an energy 
TOQ xont = a i 2 a ttA ~ IK or so - 11 a reasonable proximity 
effect can be established, then A can be on the scale 
of 1 — 10K. Fields of only a few Tesla can produce 
Zeeman energies up to ~ 100K in magnitude due to 
large spin-orbit enhancement of the <?-factor. And 
finally, t determines the bandwidth and should therefore 
be on the eV scale. These very rough estimates suggest 
that the relevant hierarchy of energies is t 3> c*i a tt > A, 
with V z on the scale of «i a tt or smaller. (Hereafter we 
cease to distinguish the lattice and continuum spin-orbit 
coupling constants.) We will adhere to this rough 
guideline in our simulations, and henceforth work in 
units where t = 1. Note, however, that as A and V z 
decrease toward zero, then finite-size effects caused by 
long coherence lengths make simulating the topological 
phase increasingly difficult. 



and ipk 



Lck (L is the system size). In real space, the 
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III. INTERACTION EFFECTS 



DMRG study of the phase diagram with 
interactions I 



As a first step towards incorporating interaction ef- 
fects, we will use DMRG to assess the phase space oc- 
cupied by the topological phase when U > 0. More pre- 
cisely, we would like to explore how the criteria displayed 
in Eq. Q, which specify where the topological phase 
appears in the non- interacting limit, evolve with inter- 
actions. Of particular experimental relevance is the re- 
quired range of /i, which is closely related to the stability 
of the topological phase against disorder-induced chemi- 
cal potential variations which will inevitably be present 
in any semiconducting wire. To protect against this kind 
of disorder, one would ideally like to engineer the system 
to remain topological over as broad a range of y, as pos- 
sible. (There are, however, competing factors associated 
with the bulk gap discussed later in this section.) Ini- 
tially, we will address these issues when A is sufficiently 
large that the wire is rather far from the gapless Luttinger 
liquid regime. In Sees. |IIIC| and |HID| we will revisit the 
problem in the opposite regime when A, V z — > 0. 

All of our DMRG simulations are performed on a finite 
system with open boundary conditions. The topological 
superconducting state can be identified numerically in 
three complementary ways. First, since fermion parity is 
a good quantum number it is straightforward in DMRG 
to determine the minimum energies £wen/odd f° r eigen- 
states in the even and odd parity sectors. The energy 
difference 



AE = \E C 



dd 



E, 



(8) 



is finite in the ordinary superconducting state (because 
there is a unique ground state) , but vanishes in the topo- 
logical phase due to the presence of Majorana modes 71,2 
exponentially localized at the left and right ends of the 
system. Thus AE provides a sort of 'order parameter' 
that distinguishes the topological and trivial supercon- 
ducting phases. 

Second, because DMRG provides direct access to the 
many-body wavefunction we may compute the bulk en- 
tanglement spectrum. Defining the left reduced density 
matrix — Ttj^l){9\ where the trace is over all sites 
in the right half of the wire, the entanglement spectrum 
consists of the energies of the entanglement Hamiltonian 
He = — In pl- In the topological phase, this spectrum 
is expected to be two-fold degenerate.^ Physically, one 
may think of the entanglement degeneracy as a precur- 
sor of the zero-energy Majorana edge states that would 
appear if one were to cut the wire through a bulk bond!^ 

In Fig. [2] we display AE (upper panel) and the entan- 
glement spectrum at the center of a 400 site chain (lower 
panel) as a function of /x with A = 0.1, a — 0.3, V z = 0.3 
and U = 0.1. The region of parameter space over which 
there is an entanglement degeneracy exactly corresponds 
to the topological phase with AE — 0. 




FIG. 2: Even/odd sector ground state energy difference AE 
(upper panel) and entanglement spectrum at the center bond 
(lower panel) for a 400 site chain with A = 0.1, a — 0.3, 
V z — 0.3 and U = 0.1. The topological superconducting phase 
is signaled by AE = together with a two-fold degeneracy in 
the bulk entanglement spectrum. 



A third way to detect the topological phase numeri- 
cally is by extracting the wavefunctions of the Majorana 
edge modes as follows. Let |0) and 1 1) be the minimum- 
energy many-body wavefunctions for the wire in the even 
and odd parity sectors. In the topological phase these 
are degenerate and can be related through the Majorana 
operators 71^ that satisfy 7^ = 1 and {7,3,7;,} = 25 a b 
according to 



|1) =e* e 7l |0) =ie i9 j2\0). 



(9) 



Here 8 is related to the overall phase factors for |0) and 
|1), and can always be absorbed into the definition of 
(say) |0). Let us assume the following expansion for 71,2: 



7a 



3 s=t,l 



h.c), 



(10) 



where (f>^ is the wavefunction corresponding to Majo- 

is given by the anti- 



rana mode 7 a . The component 
commutator 



{ C L 7a} = (Pif- 
Using Eqs. (|9| and ( 11 ), one can show that 



iff 



J6, 



i|4lo>- 
(i|4io) 



-ie 1 



ie 



- w 1 



lili) 

let ID 



(11) 

(12) 
(13) 



Thus measuring (l|qj0) and (0|cj s |l) numerically allows 
one to back out the precise form of the Majorana wave- 
functions when one is in the topological phase. Figure [3] 
shows an explicit DMRG calculation of the probability 

distribution py' —J2 S l^sl 2 01 tne Majorana wavefunc- 
tions, each of which are indeed localized at a single edge 
in the topological phase. 
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FIG. 3: Probability distribution p) a> = J2 B \<P' 
single-particle Majorana edge state wavefunctions <^ < - 1 ' (left 
panel) and (f>^ (right panel) for a 400 site chain with A = 0.1, 
a — 0.3, V z — 0.3, U = 0.1 and /j, = 0. For these parameter 
values the system is in the topological phase (see Fig.|2|, thus 
(f)' 1 ' 2 ' are non-zero only near a single edge of the system and 
decay exponentially into the bulk. 



A cautionary remark is in order. The procedure out- 
lined in the preceding paragraph relied on our assumption 
in Eq. ( 10 1 that the Majorana operators are purely linear 

This assump- 



in the lattice fermion operators Cj S and 
tion of course holds in the non-interacting limit, though 
there is no obvious reason why higher-order terms (in- 
volving, e.g., three-fermion components) should be for- 
bidden in an interacting system. One can in fact write 
down simple exactly solvable interacting models where 
the Majorana operators involve only three-fermion terms. 
Nevertheless, in all our simulations, including those with 
very strong interactions, we empirically find that such 
corrections, if present at all in our model, are exceedingly 
weak. This can be deduced by computing (j>^ and (j>^ 
as outlined above and then computing their normaliza- 
tion. In all cases we examined the normalizations deviate 
from unity by less than 1%, strongly suggesting that a 
single-body decomposition of the Majorana operators is 
indeed adequate. 

Using these methods, one can employ DMRG to deter- 
mine the parameter range in which the topological phase 
exists, for arbitrary-strength interactions. Figure [4] il- 
lustrates the phase diagram for a 400-site chain with 
A = 0.1 and a — 0.3 as a function of /i and V z , with 
three different interaction strengths. The finite-C/ phase 
boundaries were obtained by sweeping V z at fixed /i, and 
plotting the value at which AE first vanishes. Note that 
the left phase boundary changes very little with interac- 
tions since the electron density is extremely low in that 
region of the phase diagram. The minimum of the topo- 
logical phase boundary additionally shifts to finite fi; this 
property arises because the U repulsion adds a charg- 
ing energy and thus effectively renormalizes the chemical 
potential. Two more important trends are also evident 
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FIG. 4: Phase diagram showing the broadening of the topo- 
logical superconductor phase (shaded regions) with increasing 
interaction strength. The data was obtained for a 400 site 
chain with A = 0.1 and a — 0.3. 
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FIG. 5: Broadening of the topological superconductor phase 
at fixed Zeeman field of V z = 0.3 with increasing interaction 
strength. The data corresponds to a 400 site chain with A = 
0.1 and a — 0.3. Error bars are smaller than the symbol sizes. 



in the figure: as U increases, 1) the minimum value of 
V z required to stabilize the phase decreases significantly, 
and 2) at larger V z the topological phase occurs over a 
much broader chemical potential window. The latter re- 
sult is further illustrated in Fig. [5] which displays the 
phase boundary as a function of U at fixed V z = 0.3. 
Both features are desirable from an experimental stand- 
point. The decrease in the required V z in principle allows 
one to apply weaker fields, thereby disturbing the prox- 
imate superconductor less, while as already mentioned 
the widening of the topological phase as a function of /i 
implies greater immunity against chemical potential fluc- 
tuations in the wire. 

The phase diagrams shown in Figs. [2] and [5] admit a 
natural interpretation if we think of repulsive interac- 
tions as enhancing the effective Zeeman splitting while 
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FIG. 6: Energy of the first and second excited states rela- 
tive to the ground state for interaction strengths (7 = and 
U — 0.1. In the topological phase the first excited state is de- 
generate with the ground state. The data was again obtained 
for a 400 site chain with a = 0.3, A = 0.1 and V* = 0.3. 



suppressing the pairing A. Electrons can avoid an energy 
cost from the repulsion U by aligning their spins. Also, 
U has the 'wrong' sign for favoring pairing. According 
to Eqs. if introducing U > effectively enhances V z 
while suppressing A, then qualitatively similar behavior 
to what we found numerically ought to occur. However, 
this interpretation also suggests that the experimental 
benefits mentioned above are not without a cost. For ex- 
ample, the more spin-polarized the electrons are, the less 
effective the s-wave pair field A is at inducing pairing. 
Thus one expects a decrease in the bulk excitation gap 
in the topological phase. 

To test this picture we measure the effect of inter- 
actions on the bulk gap and magnetization within the 
topological phase. Because of the ground state degen- 
eracy, one must compute the lowest three eigenstates 
of the interacting Hamiltonian to capture the bulk gap. 
This poses a challenge for the usual DMRG technique of 
targeting multiple states as one would have to represent 
three states using the same basis. To overcome this dif- 
ficulty, we instead compute the lowest three eigenstates 
one at a time, constraining each state to be orthogonal 
to those previously calculated by introducing an energy 
penalty for any non-zero overlap. This requires us to be 
able to store and manipulate wavefunctions from inde- 
pendent DMRG calculations which we achieve by work- 
ing in the matrix product state formalism.^ 

Figure [6] shows how AE and the bulk gap E gap depend 
on ji for [7 = and U = 0.1, keeping the other param- 
eters fixed to A = 0.1, a — 0.3, and V z = 0.3. In order 
to be sure that the gap we obtain is truly characteristic 
of the bulk, rather than arising from an edge excitation, 
we have checked that local properties of the second ex- 
cited state, such as the electron density, are significantly 
different from the ground state throughout the system. 



The bulk gap turns out to be roughly constant across the 
topological phase, changing rapidly only near the phase 
boundaries. (In an infinite system, the bulk gap can be 
set either by excitations carrying zero momentum or mo- 
mentum near the Fermi points, depending on parame- 
ters. This accounts for the abrupt change in slope near 
/i = 0.2 in the figure.) We can therefore treat the maxi- 
mum value of the gap as a measure for the robustness of 
the topological phase (to thermal fluctuations). 

Figure [7](a) shows this maximum gap value as a func- 
tion of interaction strength for various values of the spin- 
orbit coupling. Consistent with the prediction of Gan- 
gadharaiah et alW^ for the large- V z limit and the physical 
picture discussed above, we find that the bulk gap mono- 
tonically decreases as the interaction strength increases. 
The actual values observed range from about 20% to 80% 
of the pairing field A. 

Next, to check the intuition that taking U > should 
lead to greater spin alignment, we measure the magne- 
tization per particle which, like the bulk gap, is roughly 
constant as a function of \i over most of the topological 
phase. We can therefore again take its maximum value 
as an indication of the magnetization for the phase. The 
data appears in Fig. [7Fb) and clearly indicates that in- 
teractions effectively enhance the Zeeman field. 



B. Hartree-Fock analysis 

Remarkably, qualitative and even semi-quantitative 
agreement with the phase diagram calculated from 
DMRG can be obtained within a simple Hartree-Fock ap- 
proach. Of course Hartree-Fock theory is not known for 
its spectacular success in one-dimensional systems, where 
fluctuations typically negate its predictions even on a 
qualitative level. Here, however, the presence of Spin- 
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FIG. 7: Interaction dependence of (a) the maximum bulk en- 
ergy gap and (b) the maximum magnetization per particle 
within the topological phase for a 400 site chain with param- 
eters A = 0.1 and V z = 0.3. Error bars are smaller than the 
symbol sizes. 
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orbit coupling, the Zeeman field, and pairing conspire 
to eliminate all continuous symmetries for the problem, 
rendering the system rather 'stiff' to fluctuations. Conse- 
quently Hartree-Fock theory is expected to perform rea- 
sonably well here, at least when these symmetry-breaking 
fields are sufficiently strong. 

To facilitate comparison with the DMRG results above, 
we continue to work with a lattice model, though now 
with periodic boundary conditions. We will consider a re- 
stricted set of Hartree-Fock trial states which are ground 
states of the following effective non-interacting Hamilto- 



HpfF = 



E 



3 

a 
2 ( 



i&j(T y Cj + i + h.c.) 



V z c)a z Cj 



(14) 



Notice that this has the same form as H in Eq. @, 
except that now t, fl, a, V z , and A are regarded as varia- 
tional parameters. One can diagonalize H e g by going to 
momentum space and then expressing 

c fcT = ^{k)d Ak + ^{k)d Bk + ^f{k)d\_ k + ipf{k)d B _ k 

C U = V2( k ) d Ak + ^2( k ) d Bk + ^i(k)d A _ k + (pf(k)d^ g _ k . 

(15) 

Here <$ A / Bk create quasiparticle excitations with energy 

A/ B 

E A / B (k) > and (k) are the corresponding wave- 
function components. While the energies and wavefunc- 
tions can be obtained analytically from the Bogoliubov- 
de Gennes equation, here we simply note that due to 
spin-orbit coupling and are odd in k while all 
other components are even. The ground state of H c g 
is annihilated by all d A / Bk operators and can thus be 
written in terms of the vacuum |vac) of Cj s fermions as 



= J^d^fcdsfclvac). 



(16) 



(This ground state is unique in the periodic boundary 
condition geometry.) 

The variational parameters are selected by minimizing 
the expectation value of the interacting Hamiltonian H 
with respect to the trial state 



Ehf — E + Ej 
E = (i>\H \il>) 
Ei = MHM). 



(17) 
(18) 
(19) 



The Hq expectation value is of course minimized when 
the variational parameters equal their respective values 
in the original, interacting Hamiltonian. Deviations away 
from these values may become favorable, however, when 



an increase in Eq is offset by an energy gain from Ej. 
The latter may be explicitly expressed as 



E, 



(N t N l + \A\ 2 ), 



(20) 



where N s is the average fermion number with spin s and 
A is the anomalous correlator that is non-zero due to 
pairing: 



k 



(21) 
(22) 



(The exchange term — | X)fc( c lt Cfc ^ I 2 vamsnes by sym- 
metry, which can be seen formally due to the parity 
of the wavefunction components mentioned above under 
k — > —k.) Thus there are two ways in which the system 
may lower its interaction energy: by enhancing its spin 
polarization to decrease the first term in Eq. ( 20 1 and by 



pairing less strongly to diminish the second term in Eq. 

poi. 



We deduce the outcome of the competition between 
E and Ej, and hence the fate of the topological phase, 
by numerically minimizing the Hartree-Fock energy Ehf 
over a range of V z and [i. Since the overall scale of the 
parameters in H c h does not affect Ehf, we set t = t = 1 
at the outset. We find also that the optimal a is always 
very close to a, the difference being likely due to nu- 
merical error. The remaining parameters are, however, 
significantly renormalized by interactions — fl because of 
a charging energy associated with U, and A, V z for the 
reasons described above. With the optimal parameters 
in hand, Hartree-Fock theory predicts that the system 

is topological when V z > A and —\/V z — A 2 < fi < 



V 2 — A 2 , which constitutes a straightforward general- 
ization of Eqs. Q. 

Figure [8] illustrates our results when A = 0.1 and 
a = 0.3; the data in parts (a) and (b) were obtained 
with U — 0.5, while parts (c) and (d) correspond to 
[7=1. Dashed lines indicate the topological phase 
boundary for the non-interacting system, square symbols 
denote DMRG data points from Fig. |4j and thick solid 
lines represent the Hartree-Fock phase boundary. The 
shading in Figs.pSJa) and (c) shows the optimal V z , which 
is enhanced throughout compared to the bare V z , most 
strongly inside of the topological phase and for larger 
values of /i. Similarly, the shading in Figs. [8jb) and (d) 
indicates the optimal A normalized by the bare A. This 
is suppressed throughout, though most weakly when the 
Zeeman splitting is large. The 'dual' role of the Zeeman 
field underpins this trend: V z both magnetizes the system 
and simultaneously suppresses pairing. Hence at large V z 
the anomalous correlator in Ej above is already signif- 
icantly suppressed, even when A = A. The renormal- 
ization of the Zeeman and pairing fields clearly underlies 
the expansion of the topological phase evident in Fig. 



8 



(a) 




(b) 




-0.4-0.2 0.2 0.4 

(d) 




-0.4 ,, 0.4 0.8 



FIG. 8: Hartree-Fock phase diagrams with A = 0.1, a = 
0.3, and interaction strengths U = 0.5 [parts (a) and (b)] 
and U = 1.0 [parts (c) and (d)]. The dashed lines repre- 
sent the phase boundary separating the topological and non- 
topological phases when (7 = 0, while the thick solid lines 
represent the Hartree-Fock phase boundary for an interacting 
wire. Square symbols denote DMRG data points from Fig. 
[4] The effective Zeeman field V z — which interactions enhance 
compared to V z — is shown by the shading in (a) and (c). The 
effective pairing field A normalized by A is represented by 
the shading in (b) and (d); this is clearly suppressed by inter- 
actions. 

[8j lending further support to the interpretation provided 
above for our DMRG results. 

The Hartree-Fock and DMRG phase boundaries agree 
remarkably well in the U — 0.5 case, with Hartree-Fock 
theory only marginally overestimating the extent of the 
topological phase. Poorer agreement occurs, however, 
when {7=1. Aside from the greater quantitative over- 
estimate of the width of the topological phase at finite 
V z , Hartree-Fock theory fails in a more serious way as 
V z — » 0. Namely, as Figs. 8jc) and (d) indicate, the 
topological phase is predicted to exist over a finite range 
of fj, even at V z =0. This can only arise if a finite Zee- 
man field is spontaneously generated at V z = 0, which is 
forbidden by the Mermin- Wagner theorem in the model 
under consideration (V z breaks both time-reversal sym- 
metry and a continuous U(l) spin symmetry that would 
otherwise be present in H ; see Sec. |III E| for a more de- 
tailed discussion). 

Thus while Hartree-Fock theory provides a simple, in- 
tuitive picture for the influence of interactions on the 
topological phase, its regime of validity is, not unexpect- 
edly, limited. This approach's inevitable shortcomings 
surface when interactions become too strong (compared 



to the gap in the non-interacting limit). There the stiff- 
ness of the system against fluctuations — which is what 
we argued might allow Hartree-Fock theory to work in 
the first place — is overwhelmed by the interactions. This 
limit is best attacked instead by first considering an inter- 
acting, spin-orbit coupled wire (which forms a Luttinger 
liquid), and then including the effects of A and V z . One 
can pursue this route using bosonization, which was al- 
ready done in Ref. iTTlin the large-Zeeman case where the 
problem could be mapped onto that of spinless fermions. 
In the following subsection we will consider the opposite 
limit where both A and V z are small. 



C. Bosonization analysis at small A, V z 

In the limit of small A and V z the natural starting 
point is a critical system of interacting electrons subject 
to spin-orbit coupling. To describe this regime analyti- 
cally we will focus now on the continuum Hamiltonian H 
in Eq. ([I]) rather than the discrete lattice model studied in 
the previous two subsections. The single-particle energy 
spectrum in this case is given by two horizontally shifted 
parabolas crossing at k — 0, as shown by the dashed lines 
in Fig. [I] Recall from the discussions leading to Eqs. Q 
that with A = V z = the system is closest to realizing 
the topological phase when [i = 0. In the non-interacting 
limit, perturbing this point with couplings A and V z such 
that V z > A then drives the wire into a topological su- 
perconductor. We would like to now understand how an 
interacting wire with fj, enters the topological phase 
upon adding weak pairing and Zeeman fields A, V z . 

We find it convenient for the following to first unitar- 
ily rotate the Hamiltonian in Eq. ([IJ by n/2 about the 
a x axis so that a y — > —a z and a z — > a v . (This is im- 
plemented by the unitary operator U = exp[— i7ra x /4]). 
In this basis spin- up electrons with a z = +1 live on 
the right-shifted parabola and exhibit Fermi points at 
kpi = 2ma and kp2 — 0, while spin-down electrons be- 
long to the left-shifted parabola with Fermi points at 
momenta —kp\ and kp2- Low-energy excitations near 
these Fermi points can be captured by decomposing the 
fermion operators in terms of slowly varying left- and 
right-moving fields ipp/Rs through 

= M + Me'"' 1 *- ( 23 ) 

The non-interacting part of the effective low-energy 
Hamiltonian with V z = A = fi = is then given by 

^Rs(^ iv Fdx)^R.s + ■tp{ s (iVFd x )-lpLs , 

while the non-oscillatory parts of the perturbations we 
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are interested in read 

J X 

Ma = / M^m^Li + i>L^Rx + h.c. 

J x 



(25) 
(26) 
(27) 



We will incorporate electron-electron interactions us- 
ing bosonization 47 , expressing the low-energy fermion 
fields in terms of conjugate bosonic fields (p s , 9 a as 

i>Rs = 7 — , i>Ls = i=== • (28) 



\/2 



no 



We then form symmetric (charge) and antisymmetric 
(spin) combinations of the bosonic variables via 



(<p t + n )/V2, <p„ = (<Pt-<p i )/V2, (29) 



and similarly for the S fields. The interacting Hamil- 
tonian with V z = A = fj, = 0is represented by a sum 
Ti = Hp+lia for the charge and spin sectors. The charge 
sector is harmonic, 



it - Vp 



K p {d x e p f + -L(a^ p ) 2 



(30) 



and is parameterized by a repulsive charge Luttinger pa- 
rameter K p < 1. The spin-sector Hamiltonian contains 
in addition to harmonic pieces a marginally irrelevant 
backscattering term Ht,si 



Tier 



-gbs 



K a (d x e a ) 

3r - 3l- 



1 

Kn 



(31) 



While a bosonized expression of "Hbs is available, we find 
it convenient to formulate this interaction in terms of spin 
currents J^/ L — ^Wr/l^^R/l describing spin-density 
fluctuations at the right/left Fermi points. In this rep- 
resentation the harmonic part of 1-L a describes the spin 
sector for non-interacting electrons, with all interaction 
terms for this sector encoded in Hbs- For weak interac- 
tions the backscattering amplitude is given by gbs ~ 211 . 
Also, for both sectors the product v v K v = vp remains 
unrenormalized . 

A brief digression is in order regarding the spin 
backscattering term in Eq. (31 1 and the charge Luttinger 
parameter K p . Naively, one might expect that the SU(2)- 
invariant form of Hhs above is not appropriate for our 
system because of the presence of spin-orbit coupling. 
However, this is not the case for the specific models we 
have been considering, which in fact exhibit a 'hidden' 
SU(2) symmetry that guarantees the legitimacy of Eq. 



(31 1 . This symmetry is most simply exposed in the con- 
tinuum Hamiltonian in Eq. (JlJ when V z = A = by 
defining new operators ip — exp[—imaa y x]'ip'. When 
the Hamiltonian is expressed in terms of ?/>', the spin- 
orbit term is absent while the interaction term retains ex- 
actly the same form — i.e., the transformed Hamiltonian 
exhibits SU(2) symmetry. A similar conclusion applies 
for the lattice model in Eq. (J5|, which can be mapped 
onto the pure Hubbard model with renormalized hopping 
t' = \/t 2 + a 2 upon performing a suitable unitary trans- 
formation to eliminate the spin-orbit term.^ It also fol- 
lows that, just as in the pure Hubbard model, the charge 
Luttinger parameter appropriate for our DMRG simula- 
tions satisfies 1 * 9 1/2<K P <1, even when a =^ 0. 

In bosonized language the perturbations in Eqs 
through (27 1 read 



(25) 



rLy t — 
Ha = 



T~Lp — 



V z 



sin[\/27r(^ p + 9 a )], 



r 2A 

/ — cos[v27T0J cos[v27ry CT ], 
J x ira 



/' / \/ - d x tp p . 



(32) 
(33) 
(34) 



Observe that Hy z and Ha are expressed in terms of dual 
fields relative to one another: ip p vs. 9 p and 6 a vs. ip a 
correspondingly. The scaling dimensions of the perturb- 
ing operators are 

1 \ 1 



D 



D„ 



Kn 



1 



(1 + Kp 



(35) 



1 



1 

K~, 



(36) 
(37) 



Under renormalization the corresponding couplings grow 
according to 



dl 
dA 
H 
dfx 



(3 - K p ) 



1 9bs 
4 2ttv f 



V, 



ii-Kp 1 ) 



42irv F 



Mi 



(38) 



where £ is a logarithmic rescaling factor. The contribu- 
tions to the flow equations containing the backscatter- 
ing coupling <7bs are most conveniently derived using the 
original fermionic representations in Eqs. (25), (26) and 
(31 1. One can perturbatively expand the action in these 



terms and then apply Wick's theorem to fuse fermions at 
nearby space-time points. Key technical details of this 
procedure can be found, for example, in Ref. 50, Alter- 
natively, a standard one-loop momentum-shell renormal- 
ization group analysis leads to the same results. It is 
worth noting that in our calculation the renormalization 
group flow of <7bs is standard, 

dgbs 



d( 



2-KVl 



(39) 
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so that g^ s (£) / (2ttvf) = go/(2iTVp+go£). Here 50 denotes 
initial value of the backscattering interaction. Note also 
that n, which is expressed in terms of charge fields only 
in Eq. (34), is completely unaffected by Hbs- 



by 



The rcnormalization group equations ( 38 1 are solved 



V z (£) 
A(£) 



14 exp 

A exp 

fie*, 



(3-K p 



got V/ 4 



2-KVl 



-3/4 



(40) 



where V z , A, ji on the right-hand-side stand for initial 
values of these couplings. Thus V z and 11 constitute rele- 
vant perturbations for any repulsive interaction strength, 
while A becomes irrelevant when K p < 1/3. In what fol- 
lows we will assume 1/3 < K p < 1 so that all three cou- 
plings are relevant. (For K p < 1/3 interactions are too 
strong for superconductivity — topological or otherwise — 
to take root.) 

Let us first explore the most important case of fi = 0, 
where the competition between the Zeeman and super- 
conducting perturbations is most apparent. Since both 
are relevant in the parameter range of interest, the fate 
of the system depends on which of the two renormalized 
parameters first flows to strong coupling (i.e., values of 
order the Fermi velocity vp). If the initial conditions 
are such that V z dominates, then the topological super- 
conducting state emerges, whereas if A dominates an or- 
dinary superconducting state appears. Comparing the 
renormalization group scales £y z and £a at which the 
renormalized parameters reach strong coupling, we find 
that the topological phase arises provided 



V z > cA p with j3 



3-K p + g /(4irv F ) 



3-K 



3# /(47™_f) 



(41) 



for some constant c. Observe that for repulsive interac- 
tions the exponent /? > 1, resulting in a convex phase 
boundary in the A — V z plane. Consequently, interac- 
tions allow one to access the topological phase at weaker 
Zeeman fields relative to the non-interacting case, which 
is qualitatively consistent with what we found earlier in 
our DMRG and Hartree-Fock studies at finite A. Note 



that in deriving Eq. (41) we assumed go <C vp; in this 



regime backscattering-induced corrections to V z and A 
change little in comparison with the exponential terms 
in Eqs. (40). It is interesting that despite the marginal 



irrelevance of the backscattering, at least in the weak 
coupling regime, where 1 — K p ~ U/(nvp), the spin sec- 
tor contribution to the exponent /3 is in fact comparable 
to that of the charge sector. 

Similar considerations for the competition between V z 
and [L imply that at weak A the topological phase appears 
when 

V z > c'l/if with a (3 - K p + g /(Aixv F ))/2, (42) 



where c' is a constant. Since j3' > 1, the chemical poten- 
tial window over which this state occurs broadens, which 
is again in line with our earlier numerical results. 

As an aside, we briefly contrast our results with 
the predictions of Gangadharaiah et alW^ for the large- 
Zeeman limit, where only two rather than four fields 
are needed to capture the low-energy physics. We can 
access their limit from our approach by first observing 
that at large V z the combination (f> p + 6„ is pinned by 
Ti.\r z 1 reducing the number of gapless modes by a factor 
of two. The remaining gapless mode is described by fields 
9 = (9 p — (/)&)/ V2 and 4> = (<f>p — a )/v2. Integrating out 
massive fluctuations and rewriting % in terms of 9, <j>, one 
obtains 



K(d x 9) 2 + -U<9^) 2 

A 



(43) 



with 



K 



Kl + (v a K a )/(v p K p ) + (v p /K p )(Kjv a ) + K p 2 



y/(l + KZ)(l + K p *) 



(44) 



(For simplicity, we neglected the marginally-irrelevant 
backscattering term here.) At the fixed point with only 
one gapless mode, the pairing term becomes irr elev ant 
for K < 1/2 1 S3 it i s tempting to try use Eq. (44) to 
reconcile this result with our finding that when V z — > 
and K a = 1 pairing becomes irrelevant when K p < 1/3. 
Such a comparison is difficult, however, since the Lut- 
tinger parameters K p and K a will flow away from their 
V z = values as the Zeeman field renormalizes to strong 
coupling. We simply remark that in general there is no 
reason to expect the critical interaction strength above 
which pairing becomes irrelevant to be the same at weak 
and strong V z . This opens the intriguing possibility that 
there exists a range of interactions where the topological 
phase can be accessed at weak V z , but beyond a critical 
value is replaced by a gapless state. Such a scenario is 
very different from the case of non-interacting electrons 
where the topological phase can in principle appear at 
arbitrarily large V z , and would be interesting to explore 
in future numerical work. 



D. DMRG study of the phase diagram with 
interactions II 



To verify our perturbative bosonization predictions 
and test their range of applicability we will now quantita- 
tively compare these results with DMRG data. Doing so 
requires us to extend our DMRG simulations to smaller 
values of A than we analyzed previously. To minimize 
finite-size effects caused by long correlation lengths in 
this regime, our numerics for this section were performed 
on an 800-site system which is twice as large as those 
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FIG. 9: Critical V z values separating the topological phase 
(V z > V z clit ) from the conventional phase (V z < V z clit ) . The 
DMRG data was obtained by scanning V z at fixed A for an 
800 site system with a — 0.3. We have chosen /i according to 
the relation /i = (0.0832) U, which we find empirically cancels 
the charging effect arising from the on-site repulsion and keeps 
the density fixed at its U = 0, fj, = value. The solid lines 
are fits of the form c A' 3 with the exponents j3 evaluated from 
Eq. (41 1 taking go — U/2 and K p from the Bethe ansatz. 



previously simulated. The data points in Fig. [9] illustrate 
the computed phase boundary separating the topologi- 
cal and conventional phases in the A — V z plane, with 
a = 0.3 and at several values of U. Qualitative agree- 



ment with Eq. (41 ) is immediately apparent; interactions 
indeed cause the topological phase to broaden and occur 
at weaker values of V z . 

For a more detailed comparison, we fit the critical V z 
values computed with DMRG to a power law of the form 
V 2 crit (A) = c A' 3 . Treating both c and f3 as fitting param- 
eters indeed allows one to fit all of the data from Fig. [9] 
extremely well (not shown), even those at the largest A 
considered where the validity of our bosonization analysis 
is not a priori obvious. It is interesting to ask whether 
(3 can be calculated rather than treated as a fitting pa- 
rameter, since K B can be obtained exactly for the pure 
Hubbard model. 49 Such an analysis is complicated, how- 
ever, by the influence of the spin-sector backscattering 
term on the exponent /3, particularly given that all of 
the U > curves are rather far from the regime where 
this term may be treated perturbatively. 

It is instructive to first attempt to fit the data by com- 
puting f3 with go artificially set to zero in Eq. (41 1. To 



obtain K p , we first perform the unitary transformation 
discussed in the previous section to map the V z = A = 
lattice Hamiltonian onto the pure Hubbard model with 
hopping parameter t' — \Jt 2 + a 2 . One can then numer- 
ically solve the Bethe ansatz equations appearing, for 
example, in Ref. [51] and determine the charge-sector ve- 
locity and Luttinger parameter by calculating both the 
compressibility and the response of the system to a small 
flux! 47 ! 49 ] [Note that an interaction of strength U in our 



DMRG simulations corresponds to 2U in the usual Hub- 
bard model literature because of the factor of 1/2 multi- 
plying the hopping term in Eq. Mfy .] When the exponent 
/3 is determined in this fashion, we find that the fit to our 
DMRG data is in fact rather poor. This strongly sug- 
gests that the spin backscattering term — despite being 
marginally irrelevant — in fact plays an important quan- 
titative role, as argued in the previous section. 

It is unfortunately difficult to obtain a reliable quanti- 
tative estimate for the initial backscattering coupling go 
outside of the perturbative regime (which does not apply 
here). Thus we will instead assume an ansatz go = bll, 
where b is a fitting parameter independent of U. Note 
that this is far more restrictive than treating j3 as a sepa- 
rate fitting parameter for each value of U. By computing 
/3 from Eq. (41 1 and treating b and c as a fitting param- 
eters, we obtain the fits shown by solid lines in Fig. [9] 
These curves correspond to b = 1/2 and values of c and 
(3 listed in the figure. The agreement between field the- 
oretic predictions and DMRG is rather remarkable, and 
holds for all A and U shown. This agreement is particu- 
larly encouraging since the two methods compared here 
are complementary in the sense that DMRG operates op- 
timally with short correlation lengths while bosonization 
works perturbatively in small A and V z where correla- 
tion lengths are long. By applying both methods one 
can seamlessly connect these two very different limits. 



E. Topological superconductivity at zero magnetic 

field 



The results from the preceding subsections collectively 
paint the following general picture: local Coulomb re- 
pulsion can significantly broaden the topological super- 
conducting state in phase space, allowing the topological 
regime to be accessed at smaller Zeeman fields and over 
a wider chemical potential window. Given the trend ex- 
posed by our DMRG data in Fig. [4j it is natural to ask 
whether, at still larger interaction strength U, the topo- 
logical phase might extend all the way down to V z = 0. 
We saw in Sec. |IIIB| that Hartree-Fock theory in fact 
makes such a prediction, but unfortunately dismissed this 
as an unphysical artifact on very general grounds. 

Let us discuss the underlying reason in greater detail 
here. Recall that realizing the topological phase requires 
the wire to enter a 'spinless' regime where there exists 
only a single set of Fermi points. Due to Kramer's theo- 
rem, this is fundamentally impossible in a time-reversal- 
invariant system. One can understand this in the con- 
text of the band structure in Fig. [I] by observing that 
the k = crossing of the dashed lines (which correspond 
to V z = A = 0) is protected by time-reversal symme- 
try. The applied Zeeman field serves the sole purpose 
of breaking this symmetry, but as an inessential byprod- 
uct simultaneously breaks a second, continuous symme- 
try exhibited by the Hamiltonians we have been con- 
sidering up to this point. Namely, when V z = the 
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FIG. 10: Phase diagram of a wire with spin-orbit coupling 
a = 0.3, proximity- induced pairing A = 0.1, and Ising-like 
interactions given by Eq. ( 47 1 . Because of the first-order na- 



ture of the phase transitions, reliable DMRG data could be 
obtained for systems as small as 100 sites, though larger sys- 
tems were also used to check results. Error bars are repre- 
sented by the symbol sizes. By using the same 'order pa- 
rameters' discussed in Sec. |III A] we observe a spontaneously 
generated topological superconducting phase over a large re- 
gion of the phase diagram without an applied Zeeman field. 



Hamiltonian H in Eq. Q (like its continuum version 
H) is invariant under global U(l) spin rotations about 
the y-axis, which is lifted when V z ^ 0. In principle, 
time-reversal symmetry — which is discrete — can be bro- 
ken spontaneously, though this U(l) symmetry cannot. 



We now pose the following question: how generic is 
the continuous U(l) symmetry which prevents the spon- 
taneous generation of a Zeeman field? For the strongly 
spin-orbit coupled wires we have implicitly been consid- 
ering throughout, this U(l) is in fact quite far from being 
a microscopic symmetry. Its presence merely reflects our 
inclusion of substrate-induced Rashba coupling only, and 
our neglect of Dresselhaus spin-orbit terms that are in- 
trinsic to the wire. A more realistic model which accounts 
for both types of spin-orbit interactions will generically 
exhibit, at most, time-reversal and the discrete space- 
group symmetries exhibited by the particular wire under 
consideration. The barrier to generating a Zeeman field 
spontaneously through strong interactions is then lifted, 
in principle allowing the topological phase to be accessed 
without the explicit application of a magnetic field. 



As proof of this principle, we use DMRG to study a 
lattice model where the U(l) symmetry is broken not 
by additional spin-orbit terms, but rather through the 




FIG. 11: Cut through the phase diagram of Fig. [lOjat fixed 
(i — 0.3. Once J increases to J ~ 0.693 the magnetization 
undergoes a strong first-order jump to a finite value and the 
system enters the topological phase signaled by AE — 0. At 
the upper boundary J ~ 0.872 the magnetization and site 
occupancy both saturate to a value of nearly 1/2 and the 
system undergoes a transition to a conventional time-reversal 
breaking superconducting phase. 



interactions themselves: 



H' 



K+H'x (45) 



H'j = -jJ2(cy Cj )(cl +1 a z c j+1 ), 



(46) 
(47) 



with J > 0. Here H' describes a Rashba-coupled wire 
with proximity-induced pairing, but without a Zeeman 
field. The interaction term Hj favors ferromagnetically 
aligning the spins either along or against the z direction, 
thereby lifting the U(l) spin symmetry exhibited by Hq. 
Note that the entire Hamiltonian is time-reversal sym- 
metric and possesses no continuous symmetries. 

Before turning to DMRG to determine this model's 
phase diagram, it is useful to consider some simple lim- 
its with a, A finite. In the non-interacting case J = 0, 
the absence of a Zeeman field implies that the Hamilto- 
nian can only capture an ordinary, non-topological su- 
perconducting phase. Since this phase is protected by a 
finite energy gap it must persist up to some finite value 
of J. In the opposite limit where J — ¥ oo, the model re- 
duces to a ferromagnetic Ising chain which spontaneously 
breaks time-reversal symmetry by singly occupying ev- 
ery site with perfectly spin-polarized electrons. For J 
large (compared to all energy scales) but finite, spin- 
orbit coupling leads to imperfect spin-polarization and 
the pairing field then renders the system weakly super- 
conducting. Though time-reversal symmetry breaking 
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still survives here, this superconducting state is neverthe- 
less non-topological because there are no partially filled 
bands. 

The spontaneous onset of topological superconductiv- 
ity requires both time-reversal symmetry breaking and a 
particular window of non-interacting electron densities. 
While it is far from obvious that these conditions are 
ever met in this model, our DMRG results presented in 



Fig. [10] reveal that topological superconductivity exists 
over a significant region of the phase diagram. In Fig. [11] 
we show the characteristic behavior of the system along a 
cut through the phase diagram at fixed fi = 0.3. As J in- 
creases from zero, the system undergoes a strongly first- 
order transition beyond which the electron spins align 
ferromagnetically; AE simultaneously drops to zero sig- 
naling the onset of the topological phase. At the upper 
boundary, the magnetization saturates to a value that is 
close to full spin polarization and the system becomes 
nearly half filled. 



Though we do not expect the Ising interaction Eq. ( 47 ) 



to arise in an actual experimental wire, we take the ex- 
istence of topological order in this model as motivation 
for future investigation of more realistic models that may 
spontanously enter a topological phase. For instance, one 
might consider additional spin-orbit terms that break the 
problematic U(l) spin symmetry together with strong 
Coulomb interactions. We leave such considerations for 
a future study. 



IV. DISCUSSION 

Interacting wires proximate to a bulk s-wave supercon- 
ductor provide an experimentally relevant research area 
in which Coulomb repulsion, magnetism, and Cooper 
pairing coexist and compete in interesting ways. More- 
over, because the low-energy degrees of freedom reside in 
one dimension, the problem is amenable to powerful nu- 
merical and analytical approaches including DMRG and 
bosonization, and interestingly can even be attacked rea- 
sonably well using far cruder Hartree-Fock theory when 
the proximity effect and applied magnetic fields are suf- 
ficiently strong. 

Here we employed all three of these methods, focusing 
on how on-site repulsion affects the onset of topological 
superconductivity in the wire. A rather intuitive pic- 
ture emerged from these studies — interactions suppress 
pairing and hence the bulk gap protecting the topolog- 
ical phase. As a corollary, the localization length of 
the Majorana modes in this state increases with inter- 
actions. While these features are undesirable from an 
experimental standpoint, there is an important silver lin- 
ing. Namely, interactions additionally enhance the mag- 
netization in the wire, thereby broadening significantly 
the chemical potential window over which the topologi- 
cal phase emerges. This implies a greater robustness of 
the topological phase against chemical potential fluctua- 



tions. In our view this boon more than compensates for 
the decrease in the bulk gap, since it is likely that — at 
least initially — disorder in the wire will present the main 
obstacle to be overcome as opposed to low-temperature 
requirements. We further proposed that in wires with 
appreciable Rashba and Dresselhaus coupling, strong re- 
pulsive interactions may eliminate the need for applied 
magnetic fields altogether. As proof of principle, we 
demonstrated using DMRG the emergence of topological 
superconductivity in a strongly interacting, time-reversal 
invariant toy model. 

Many interesting future directions remain to be ex- 
plored. For instance, including longer-range repulsion in 
DMRG simulations will allow one to access the regime 
in which the charge Luttinger paramete r K p d rops below 
1/2. As discussed in Ref . ST] and Sec. IIIC here inter- 



actions can have qualitatively new effects — notably, the 
complete destruction of gapped superconducting states — 
that are inaccessible with only on-site repulsion. Study- 
ing the crossover between the weak- and strong-Zeeman 
limits where the condition for destroying superconduc- 
tivity can be rather different would be particularly inter- 
esting. Extending such simulations to incorporate both 
Rashba and Dresselhaus coupling is another worthwhile 
direction to pursue, and may allow one to find more real- 
istic models where strong interactions lead to the forma- 
tion of Majorana fermions without the presence of mag- 
netic fields. Multi-channel wire s, wh ich have received 
considerable attention recentlyj^ " 32 * 52 ^ ^ can also be effi- 
ciently studied with DMRG by considering multi-leg lad- 
ders. Finally, a conceptually very interesting problem rel- 
evant for cold-atom realizations of one-dimensional topo- 
logical superconductivity 2 ^) is the prospect of seeing rem- 
nants of Majorana physics in an attractively interacting 
wire, without a proximate superconductor. This poses 
quite a subtle problem because here pairing correlations 
will die off as a power law, but should nevertheless be 
tractable using DMRG. 
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